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Abstract 

Background and methods: A longitudinal Anopheles gambiae si. insecticide-resistance monitoring programme was 
established in four sentinel sites in Burkina Faso. For three years, between 2008 and 2010, WHO diagnostic dose 
assays were used to measure the prevalence of resistance to all the major classes of insecticides at the beginning 
and end of the malaria transmission season. Species identification and genotyping for target site mutations was also 
performed and the sporozoite rate in adults determined. 

Results: At the onset of the study, resistance to DDT and pyrethroids was already prevalent in An. gambiae si. from 
the south-west of the country but mosquitoes from the two sites in central Burkina Faso were largely susceptible. 
Within three years, DDT and permethrin resistance was established in all four sites. Carbamate and 
organophosphate resistance remains relatively rare and largely confined to the south-western areas although a 
small number of bendiocarb survivors were found in all sites by the final round of monitoring. The ace-f target site 
resistance allele was present in all localities and its frequency exceeded 20% in 2010 in two of the sites. The 
frequency of the 1014F kdr mutation increased throughout the three years and by 2010, the frequency of 1014F in 
all sites combined was 0.02 in Anopheles arabiensis, 0.56 in An. gambiae M form and 0.96 in An. gambiae S form. This 
frequency did not differ significantly between the sites. The 1014S kdr allele was only found in An. arabiensis but its 
frequency increased significantly throughout the study (P = 0.0003) and in 2010 the 70745 allele frequency was 0.08 
in An. arabiensis. Maximum sporozoite rates (12%) were observed in Soumousso in 2009 and the difference 
between sites is significant for each year. 

Conclusion: Pyrethroid and DDT resistance is now established in An. gambiae si. throughout Burkina Faso. Results 
from diagnostic dose assays are highly variable within and between rounds of testing, and hence it is important 
that resistance monitoring is carried out on more than one occasion before decisions on insecticide procurement 
for vector control are made. The presence of 70745 in An. gambiae si., in addition to 7074F, is not unexpected given 
the recent report of 70745 in Benin but highlights the importance of monitoring for both mutations throughout the 
continent. Future research must now focus on the impact that this resistance is having on malaria control in Burkina 
Faso. 
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Background 

Insecticide resistance monitoring should be an essential 
component of any malaria vector control programme. 
Vector control activities have been dramatically scaled up 
in recent years across much of sub-Saharan Africa with a 
concomitant reduction in malaria cases [1]. In order to 
maintain this success, insecticide resistance management 
plans must be incorporated into National Malaria Control 
Programmes (NMCP). Burkina Faso distributed approxi- 
mately eight million long-lasting insecticide-treated nets 
(LLINs) in the past 16 months and indoor residual spray- 
ing (IRS) with bendiocarb has been ongoing since 2009 in 
the south-western part of the country. Although several 
published papers have indicated the presence of resistance, 
or suspected resistance, to all the major classes of insecti- 
cides used in malaria control in some sites in Burkina Faso 
[2-6], routine insecticide resistance monitoring is not yet 
integrated in the NMCP. 

The pattern of emergence of insecticide resistance in 
Burkina Faso mirrors that in many other countries, with 
small foci of resistance expanding rapidly in the past 
decade. Resistance to pyrethroids was detected in the 
south-western region of Burkina Faso as early as 1999 
[7,8] and remained restricted to this region, which is the 
site of the majority of the cotton growing and other 
commercial agriculture, when next surveyed in 2006 [9]. 
The frequency and distribution of the target site muta- 
tion, 1014F, conferring resistance to pyrethroids and 
DDT, increased during this time period. At the beginning 
of the century, 1014F was found only in the south- 
western region and largely restricted to the Anopheles 
gambiae S form [10] but by 2006, 1014F was found in 
An. gambiae M and S form and Anopheles arabiensis 
and had spread to the central region of the country. 

The ace-l R mutation, associated with organophosphate 
and carbamate resistance, has also been detected in An. 
gambiae s.l. in the south-western region of Burkina Faso 
[3,5]. 

The WHO/TDR network on insecticide resistance in 
African malaria vectors was established in 2008 and ran 
for three years in four countries: Benin, Burkina Faso, 
Chad and Sudan [6]. One of the major objectives was to 
improve the monitoring of insecticide resistance in areas 
where malaria control programmes are being implemen- 
ted. Here, results from the three-year study in Burkina 
Faso are reported. They demonstrate that pyrethroid re- 
sistance is now present in all four sentinel sites and that 
both the 1014F and 1014S kdr alleles are firmly estab- 
lished in An. gambiae s.l. throughout the country. 

Methods 

Study area and insecticide bio-assays 

Mosquito larvae were collected in four localities 
throughout the country (Figure 1) according to the agro- 



ecological areas and the intensity of insecticide use in 
agriculture or for public health purpose. 

Goundry (N 12° 30, W 01° 20) is in the Sahelian zone, 
in the south of Ouagadougou, 30 km far from the capital 
city. There is a low insecticide usage for agricultural pur- 
poses and this area was not targeted by the NMCP for 
LLIN distribution before 2010. 

Koupela (N 12° 10, W 00° 21) is in the Sudan Sahelian 
zone, approximately 130 km east of the capital city, 
Ouagadougou. This village is characterized by gardening, 
rice cultivation and food crops. Extensive distribution of 
LLINs has been completed by the NMCP and other social 
partners. 

Kuinima (N 11° 09, W 04° 17) is in the Sudan Savannah 
zone, south-west Burkina Faso, approximately 4 km from 
Bobo Dioulasso. This site is characterized by food crops 
and gardening where insecticides are used to protect vege- 
table crops from insect damage. No large scale distribution 
of LLINs was carried out pre- 20 10. 

Soumousso is a rural village located (N 11° 00, W 04° 
03) in the Sudan Savannah zone, in the south-western part 
of the country. It is situated 38 km from Bobo Dioulasso 
and is characterized mainly by rice and cotton cultivation. 
Pyrethroids, organophosphates and carbamates are widely 
used to control agricultural pests. As in Kunima, this area 
was not targeted for LLIN distribution before 2010. 

Collections were made twice each year, at the end of 
July and end of October, from 2008 to 2010. Mosquitoes 
were collected as larvae and reared to adults for insecti- 
cide bio-assays as described in Ranson et al. [6]. 

Species and molecular form identification 

Genomic DNA was extracted from individual mosquitoes 
using a modification of the Livak protocol [11]. In 2008 and 
2009, a random sample of between 36 and 44 mosquitoes 
from the bio-assay controls, plus all the survivors of insecti- 
cide and an equal number of dead specimens were identi- 
fied to species and molecular forms of An. gambiae s.s. 
using the methods described by Scott et al. [12] and Favia 
et al. [13] respectively. In 2010, the protocol was altered so 
that, rather than extracting DNA from the control bio- 
assays, two test bio-assays from each of the permethrin, del- 
tamethrin and DDT bio-assays were selected at random 
from each site in each collection round and all survivors 
and dead were used for DNA extraction. Thus in 2008 and 
2009, the sample size for species identification was between 
50 and 70 per site per collection round, whereas in 2010 
this was increased to over 150. In 2010, the species identifi- 
cation was performed at the Liverpool School of Tropical 
Medicine and the SINE method was used [14]. 

Detection of kdr alleles 

Mosquitoes from the control bio-assay (or the entire 
subset from selected bio-assays for 2010) were used to 
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determine the frequency of the 1014F and 1014S alleles 
in the general population. For these assays the TaqMan 
genotyping method was used [15]. To determine the as- 
sociation between kdr genotype and the probability of 
surviving the WHO bio-assay, the surviving mosquitoes 
from the permethrin, deltamethrin and DDT bio-assays, 
plus an equal number of dead mosquitoes were geno- 
typed using the allele specific PCR method [16,17]. 

Ace-1 R mutation detection 

Control mosquitoes from 2008 (n = 107) were genotyped 
for the ace-l R mutations using the protocol of Weill et 
al [18]. In 2010, the mosquitoes used in fenitrothion 
bio-assays (all survivors and dead from a random selec- 
tion of bio-assay experiments) were used for ace-l R 
genotyping (n = 661). 

Sporozoite detection 

Adult mosquitoes were collected from the four sentinel 
sites using a mixture of pit traps, exit traps and indoor 
pyrethroid spray catches. The heads and thorax of the fe- 
male mosquitoes were tested by ELISA for the presence 
of Plasmodium falciparum circumsporozoite protein 
(CSP) [19]. 

Statistical analysis 

Bio-assay data were considered for each insecticide, local- 
ity, collection round and year. Mortality was calculated as 



the percentage of individuals that died within 24 h of ex- 
posure. Bio-assay outcomes were assessed according to 
WHO [20]. Those with an overall mortality > =98% were 
considered susceptible, those with mortality <98% but 
>80% were considered potentially resistant, and, those 
with mortality <80% were strongly suspected to be resist- 
ant. The 95% confidence limits of the mortality were calcu- 
lated using Microsoft Excel and bio-assay mortalities from 
the first (2008) and third (2010) years were compared 
using Z proportion tests. For this analysis, data from both 
rounds of testing in each year were pooled. 

Temporal and spatial variations in members of An. 
gambiae species complex were investigated for each site 
using Fisher Exact test. 

The chances of surviving DDT or permethrin diagnos- 
tic doses for each species/molecular form were analysed 
using a two-tailed Fisher exact test to compare the prob- 
ability of survival in alive and dead mosquitoes between 
pairs of species. 

The change in the 1014F, 1014S and ace-l R frequencies 
were compared between sites and between years using 
Fisher exact test with significance level set at 0.05. The 
correlation between 1014F genotype and survival to in- 
secticide was determined using x 2 tests. 

The An. gambiae sporozoite infection rate and 95% 
confidence limits were calculated in each site. The differ- 
ence in the infection of mosquitoes has been compared 
using Z test to compare proportions between sites and 
between years for the same site. 
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Results 

Insecticide bio-assays 

The results of the insecticide bio-assays are shown in Figure 2 
(data from 2008 have been published previously [6], but are 
included here for comparison). Using the WHO definitions 
of resistance, the two populations from the arid savannah 
zone in central Burkina Faso, Goundry and Koupela, were ei- 
ther fully susceptible or showed a very low prevalence of re- 
sistance to all five insecticides tested at the start of the study 
in 2008. Three years later, resistance to DDT and permethrin 
was established in these populations. In the Sudan savannah 
region in the south-west of the country (sites Kuinima and 
Soumousso), low mortality with DDT and permethrin were 
observed throughout the study. The deltamethrin bio-assays 
fluctuated around the resistance cut-off of 80% throughout 
the three years in all four sites but with most data sets not 
meeting the WHO definition of resistance. 

Survival after exposure to the carbamate bendiocarb also 
fluctuated around the 80% level in Soumousso and Kui- 
nima. Higher mortality rates were obtained in Goundry and 



Koupela but some mosquitoes were surviving the diagnostic 
dose in these sites at the end of the study. Exposure to the 
organophosphate, fenitrothion resulted in near 100% mor- 
tality in the savannah zone but there was some indication of 
resistance developing in the south-west of the country. 

DDT mortality in Goundry and Koupela significantly 
decreased throughout the study while it increased in Kuinima 
and remained unchanged in Kuinima and Soumousso 
(Table 1). There was no significant change in the deltamethrin 
mortality except for Koupela where the mortality decreased. 
Bendiocarb mortality decreased significantly between the be- 
ginning and end of the study for Goundry and Koupela sites, 
while fenitrothion mortality increased significantly in one site, 
Soumousso (Table 1). The difference in mortality rate between 
sites for each year and each insecticide is statistically signifi- 
cant except for deltamethrin in 2010 (Table 1). 

Pre-exposure to the synergist, piperonyl butoxide largely 
restored susceptibility to pyrethroids when tested in 2010 
(Additional file 1) with mortality rates ranging from 91.3% 
in Kuinima to 98.1% in Goundry. Full susceptibility to 
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Figure 2 Insecticide bio-assay results for Anopheles gambiae s.l from 2008-2010 in two rounds of monitoring in four sentinel sites in 
Burkina Faso (data from 2008 have been already published [6] and are included for the purpose of comparison). The percentage 
mortalities 24 hours following a one-hour exposure to the WHO diagnostic dose of insecticide (with 95 % CI) are shown. The minimum sample 
size for each test was 100 non-blood fed female mosquitoes, three to five days old. 
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Table 1 Insecticide bio-assay results for Anopheles gambiae s.l. in 2008 and 2010 in four sentinel sites in Burkina Faso 
(data from the two rounds in each year have been pooled for this analysis) 



Insecticide 2008 201 0 P value (2008 vs 201 0) 

% mortality (95 % CI) % mortality (95 % CI) 



DDT 


Goundry 


o r (~ic\ ai \ 

85 (79 - 91) 


45 (37 - 52) 


<0.001 




Koupela 


85 {/ / ~ 88) 


ao (oa r~)\ 

43 (34 - 52) 


<0.001 




Kuinima 


16 (5 - 27) 


5 I (2 I - 4UJ 


1 




Soumousso 


32 (21 - 44) 


O A f^lC A A\ 

34 (25 - 44) 


0.620 




P value (between sites) 


<0.001 


0.001 




Permethrin 


Goundry 


80 (74-85) 


68 (57 - 79) 


0.003 




Koupela 


"7Q f~71 O A\ 

/o {/I - 84) 


rr\ ( a~i ~70\ 

oU (4/ - /3) 


<0.001 




Kuinima 


4o (z/ - oyj 


D4 (4z - OOJ 


u.o/y 




Soumousso 


55 (48 - 62) 


54 (44 - 63) 


0.464 




P value (between sites) 


<0.001 


0.01 3 




Deltamethrin 


Goundry 


89 (85 - 93) 


86 (78 - 94) 


0.180 




Koupela 


92 (8/ - 96) 


~70 f~7 1 O /1\ 

/o (/ 1 - 84) 


<0.001 




Kuinima 


oo f~7~7 nc\ 
88 (// - 98) 


Qr (~jr CY7\ 

86 (/5 - 9/) 


0.287 




Soumousso 


77 (66 - 87) 


84 (76 - 91 ) 


0.967 




P value (between sites) 


<0.001 


0.056 




Bendiocarb 


Goundry 


QQ (QQ 1 OCW 
yy [yo — I UUJ 


yj [yZ — yy) 


n m n 

U.U I u 




Koupela 


99 (97 - 101) 


95 (91 - 99) 


0.027 




Kuinima 


77 (64 - 90) 


80 (73 - 87) 


0.776 




Soumousso 


69 (60 - 78) 


78 (72 - 85) 


0.989 




P value (between sites) 


<0.001 


<0.001 




Fenitrothion 


Goundry 


100 (-) 


99 (97 - 100) 


0.248 




Koupela 


100 (-) 


99 (98 - 100) 


0.499 




Kuinima 


93 (89 - 98) 


95 (92 - 98) 


0.402 




Soumousso 


87 (79 - 94) 


96 (94 - 99) 


0.001 




P value (between sites) 


<0.001 


0.050 





deltamethrin was restored for the two sites tested (Goundry 
and Koupela). 

Spatial and temporal distribution of Anopheles gambiae 
complex members 

A total of 5,541 An. gambiae complex mosquitoes from lar- 
val collection were identified to species and molecular form 
by PCR. Only control samples from each year (2,644 in 
total) were used for spatial and temporal distribution. The 
results are summarized in Figure 3 and Additional file 2. In 
2008 and 2009, An. arabiensis was the most abundant spe- 
cies in the central sites of Goundry and Koupela and An. 
gambiae s.s. was more abundant in the south-western sites 
of Kuinima and Soumousso. However, by 2010, the propor- 
tion of An. arabiensis had fallen to less than one third in all 
four sites. The number of individuals with M/S hybrid pat- 
terns detected was low throughout the three years, around 
1% in all the sites with largest confidence limits in Goundry 
(0-4%). 



Within An. gambiae s.s., the S form dominated in Kui- 
nima and Soumousso in all years, whereas the relative 
frequency of the M and S forms fluctuated more widely 
in Koupela and Kuinima. In 2008 and 2009 significant 
differences in species composition were observed be- 
tween the two rounds of testing in all sites except in 
Goundry (2008) and Koupela and Kuinima (2009) 
(P < 0.05). However, when the sample size was increased 
in 2010 (n = 326 in 2008, 306 in 2009 and 2012 in 2010), 
the only significant difference (P = 0.000) was observed 
in Koupela, where a higher proportion on M form An. 
gambiae were present early in the transmission season 
than at the end of the season (Additional file 2). 

The change in sampling strategy for determining the 
species distribution in 2010 meant it was possible to use 
the same sample set to determine the probability of surviv- 
ing the diagnostic dose for each species in each site 
(Figure 4). The M and S form mosquitoes are significantly 
more likely to survive permethrin or DDT exposure than 
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Figure 3 Species identification within Anopheles gambiae s.l. for the four sentinel sites in 2008-2010. Data are presented as proportions of 
the total for each species, An. orabiensis (A), An. gambiae M form (M), An. gambiae S form (S) and the hybrid of M and S forms (H), by year and by 
site. Sample sizes are minimum of 75 per site for 2008 and 2009 and a minimum of 309 for 2010. 



An. arabiensis (P < 0.05) except for the site of Kuinima 
where no significant difference has been observed between 
the M form and An. arabiensis for permethrin exposure 
(P =0.072 and 1 respectively for M and S compared to An. 
arabiensis). The S form exhibited higher survival rate to 
DDT compared to the M form in Goundry and Kuinima 
and only in Soumousso for permethrin. 

The probability of surviving permethrin exposure for 
each of the three species/forms did not differ significantly 
between sites (S form P = 0.069, M form P = 0.337, An. 



arabiensis P = 0.053). Similarly no difference is observed 
between sites for probability of surviving DDT exposure 
(P = 0.125 for S-form, P = 0.921 for the M-form, P = 0.940). 



Frequency and distribution of 1014F and 1014S kdr alleles 
in Burkina Faso 

The frequency of the 1014F and 1014S alleles in the con- 
trol' bio-assay populations is shown in Table 2 and allele 
frequencies for 2010 are shown in Figure 5. 
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Figure 4 The percentage survival of Anopheles gambiae s.l. species/form 24 hours following exposure to DDT or permethrin. Data for 
2010 and are stratified by species for each sentinel site. Statistical comparison between molecular forms and localities were performed using Z-test 
with significant differences (P < 0.05) showed in the figure (*). 
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Table 2 The frequency of 1014F and 70745 kdr alleles in Anopheles gambiae M form (M), S form (S), M and S hybrids 
(H) and Anopheles arabiensis (B) in four sentinel sites in Burkina Faso, from 2008-2010 

2008 2009 2010 



Locality Species N Freq Freq N Freq Freq N Freq Freq p value p value 

1014F WUS 1014F WUS 1014F 1014S 1014F WUS 
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Koupela 


R 
D 
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u.uu 
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A 1 


U.Uj 


n n i 
U.Uz 
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U.Uz 
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U. I Z^f 
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o 






o 






4 


0.88 




1 






M 


14 


0.36 


0.00 


15 


0.23 


0.00 


107 


0.51 


0.00 


0.056 






S 


6 


0.92 


0.00 


3 


0.83 




89 


0.92 


0.00 


0.348 




Kuinima 


B 


37 


0.04 


0.05 


28 


0.02 


0.04 


84 


0.04 


0.07 


1 


0.879 




H 


0 






1 


1.00 


0.00 


1 


1.00 




1.00 






M 


3 


0.33 


0.00 


8 


0.69 


0.00 


35 


0.59 


0.00 


0.6845 






S 


36 


0.88 


0.00 


24 


0.96 


0.00 


192 


0.95 


0.00 


0.2437 




Soumousso 


B 


16 


0.00 


0.00 


24 


0.50 


0.02 


89 


0.03 


0.07 


0.6672 


0.217 




H 


1 


1.00 


0.00 


0 






0 












M 


4 


0.50 


0.00 


9 


0.33 


0.00 


63 


0.60 


0.00 


0.05451 






S 


51 


0.86 


0.00 


24 


1.00 


0.00 


158 


0.98 




0.000 





Frequencies of 1014F and 1014S for the three years are compared for each molecular form in each site using Fisher exact test. (-) represents value of the data 
which cannot be calculated. 



In 2010, the 1014F allelic frequency for the four sites 
ranged from 0.92-0.98 in the S form, 0.51-0.6 in the M 
form and 0.01-0.04 in An. arabiensis. There was no sig- 
nificant difference in 1014 F frequency between the four 
sites (P > 0.08) (Figure 5). 



The 1014S form was detected in An. arabiensis only. It 
was found at a low frequency (Table 2), (maximum allele 
frequency = 0.12 in 2010). The frequency increased in all 
sites between 2008 and 2010 although this change was 
only significant in Koupela. 
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Relationship between kdr genotype and bio-assay data 

The presence of the 1 01 4F allele is clearly linked to DDT 
and pyrethroid (permethrin and deltamethrin) resistance 
in An. gambiae s.l (P < 0.001) (Table 3). The samples 
size for the 1014S mutation was too low for statistical 
analysis. However, of the 60 L1014S heterozygotes, only 
11 survived exposure to either pyrethroid or DDT, sug- 
gesting that a single copy of this mutation confers little 
protection against exposure to the WHO diagnostic 
doses. 

Frequency and distribution of ace-l R allele 

The ace-l R allele was only detected in S form An. gam- 
biae from Kuinima and Soumousso in 2008 at an allele 
frequency of 0.2 in both sites. By 2010, this allele was 
found in all sites and a single M form individual with the 
ace-l R allele was also detected in Goundry (Table 4). No 
ace-l R alleles were detected in An. arabiensis. The pro- 
portion of mosquitoes containing the ace-l R allele was 
higher in the bendiocarb or fenitrothion survivors than 
the dead progeny (0.03 vs 0.41, P = 2e -16). 

Sporozoite rates 

The presence of sporozoites in the head or thorax of 
adult mosquitoes was determined as part of a wider 
study to explore the association between insecticide re- 
sistance and transmission. Analysis is ongoing but the 
crude sporozoite rates are presented here to indicate the 
relative levels of malaria transmission in each site 
(Figure 1 and Table 5). The highest frequency of sporo- 
zoite positive An. gambiae s.l. was found in Soumousso 
in each year with 46 of 378 (12% CI [5.9%-15.3%]) mos- 
quitoes positive in 2009. The sporozoite infection rate 
differed significantly between sites and a significant dif- 
ference between years was observed in the south-western 
sites (Kuinima and Soumousso) (Table 5). 



Table 3 Correlation between 1014F and 1014L genotypes 
and survival to insecticides 





LL 


LF 


FF 


X2 


P value 


DDT 












Alive 


24 


36 


201 


184.1 


0.000 


Dead 


106 


14 


23 






Deltamethrin 












Alive 


6 


11 


45 


20.6 


0.000 


Dead 


131 


40 


160 






Permethrin 












Alive 


16 


14 


117 


75.5 


0.000 


Dead 


115 


41 


84 







LL, FF represents respectively the homozygotes for 1014L and 1014F and LF 
represent the heterozygotes. Data from all sites and all species tested in 2010 
were pooled for this analysis. 



Table 4 Frequency (Freq) of ace-1 R alleles in Anopheles 
gambiae M form (M), S form (S), M and S hybrid (H) and 
Anopheles arabiensis (B) in four sentinel sites in Burkina 
Faso from 2008-2010 







2008 




2010 






Locality 


Species 


N 


Freq ace- 7 


N 


Freq ace- 1 


P value 


Goundry 


B 


20 


0 


53 


0 






H 


1 


0 


4 


0 






M 


5 


0 


1 7 


0.03 


1 .000 




S 


9 


0 


67 


0.08 


0.341 


koupela 


B 


21 


0 


68 


0 






H 


0 




4 


0 






M 


7 


0 


42 


0 






S 


4 


0 


42 


0.12 


0.562 


Kuinima 


B 


34 


0.01 


64 


0 


0.347 




H 


0 




5 


0.10 


1.000 




M 


6 


0 


8 


0 


1.000 




S 


97 


0.20 


121 


0.22 


0.058 


Soumousso 


B 


30 


0.02 


42 


0 


0.417 




H 


1 


0 


2 


0 






M 


11 


0.05 


17 


0 


0.500 




S 


170 


0.18 


120 


0.21 


0.222 



Allelic frequencies (N) are compared using a Fisher exact test. (-) represents 
values the data which cannot be calculated. 



Discussion 

This data set provides a comprehensive picture of the 
changing pattern of insecticide resistance in Burkina 
Faso over a three-year period. The information is of clear 
operational importance for malaria control in Burkina 
Faso but, in addition, it is hoped that several of the con- 
clusions drawn from this data set will be of value for 
other countries and perhaps stimulate further debate 
about how resistance is measured in the field. 

Focusing firstly on the pattern emerging from the differ- 
ent sentinel sites it is clear that resistance to DDT and per- 
methrin is now firmly established in all four sites. At the 
beginning of the study, high survival rates were obtained in 
the south-western savannah zone sites whereas no resist- 
ance was detected in the central Sahelian sites. However, at 
the end of three years, resistance to these insecticides had 
clearly been established in An. gambiae s.l. in all four sites 
and hence, the distribution of resistance has been extended. 
It is important to note that nothing can be concluded about 
the level or impact of resistance from this data set. All that 
can be concluded is that a greater proportion of An. gam- 
biae s.l. are surviving the diagnostic dose but, as will be dis- 
cussed later, this can be partially attributed to a change in 
the composition of the species complex in some of the 
sites. To address the question of whether resistance is on 
the rise, data on whether the mosquitoes are becoming 
more resistant, ie, able to withstand higher concentrations 
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Table 5 Sporozoite rates in adult collections from study sites 



Locality Year 





2008 




2009 




2010 




P value 




Sporozoite rate 


N 


Sporozoite rate 


N 


Sporozoite rate 


N 




Goundry 


0.029 [0 - 0.045] 


175 


0.046 [0.008 - 0,011] 


151 






0.395 


Koupela 


0.014 [0 - 0.035] 


213 


0.038 [0.012-0.051] 


397 


0.05 [0.041 - 0.084] 


220 


0.115 


Kuinima 


0.008 [0 - 0.008] 


130 


0.043 [0 - 0.067] 


47 


0 [0 - 0] 


129 


0.038 


Soumousso 


0.060 [0 - 0.067] 


151 


0.122 [0.059 - 0.153] 


378 


0.061 [0.062 - 0.086] 


212 


0.015 


P value 


0.025 




< 0,0001 




0.020 







Proportion of Circumsporozoite Protein {CSP)-positive Anopheles gambioe s.l., total numbers of tested mosquitoes and P value of proportions compared between 
years for each study site and P values for comparison of infected mosquitoes proportions between locality for the same year. (-) represents value of the data 
which cannot be calculated 



of insecticide or longer exposure times is needed. This in- 
formation is not obtainable from diagnostic dose assays 
and dose or time response assays need to be performed. At 
present, very little information is available about the rela- 
tionship between surviving exposure to the diagnostic dose 
assays and the efficacy of interventions such as LLINs or 
IRS. Thus, while diagnostic doses are very useful for com- 
paring between sites, as resistance becomes established in 
populations, more information is needed in order to under- 
stand the epidemiological impact of this resistance. 

A second observation that may have resonance for 
other African countries is the differing response of mem- 
bers of the An. gambiae complex to pyrethroids or DDT. 
In all four sites the S form had the best chance of surviv- 
ing DDT or permethrin exposure, followed by the M 
form and lastly An. arabiensis, which was the most sus- 
ceptible. In fact, by the end of the study, the chance of 
surviving DDT or permethrin exposure was equal for all 
sites when broken down into species. Hence, under- 
standing species dynamics is essential for predicting the 
spread of resistance. As an example, the increase in per- 
methrin and DDT resistance in Goundry and Koupela 
during the course of the study is at least partially 
explained by the decreased proportion of An. arabiensis 
in these sites in 2010 (Figure 3). An important implica- 
tion of this observation is that, in sites where multiple 
members of the An. gambiae complex are present, resist- 
ance must be monitored at different periods throughout 
the transmission season to capture this variation. 

Riehle et al [21] reported a very high frequency of 
individuals showing an M/S hybrid pattern in two of the 
sites in the current study (Goundry and Koupela) in 

2008 and 2009, and proposed that these belonged to a 
new subspecies, named Goundry. Interestingly, although 
an unexpectedly high number of hybrids were detected 
in Goundry in 2008 (4.8%), no M/S hybrid patterns were 
observed in the 387 samples genotyped from Goundry in 

2009 and 2010. The overall frequency of M/S hybrids in 
Koupela was approximately 1% for the three years, the 



reason for this discrepancy warrants further 
investigation. 

It is unclear why the prevalence of resistance should differ 
so much between forms considering that all adults bio- 
assayed were reared from sympatric larvae and hence could 
be expected to be under the same selective pressure. Of 
course, insecticide resistance is a dynamic process and it is 
conceivable that within a matter of years resistance will be 
prevalent in all members of the species complex but the dif- 
fering rates at which this is emerging in different members 
of the An. gambiae complex remains unexplained. 

Equally intriguing is the failure for 1014F to become 
fixed in the S form of An. gambiae despite fluctuating 
around the 0.9 allele frequency for three years. Is a fit- 
ness cost counteracting the selection pressure? 

The establishment of the 1014S allele in An. arabiensis 
in West Africa has already been noted in Benin during 
the TDR project [22]. The first report of this allele in the 
An. gambiae complex was in An. gambiae s.s. in Kenya 
where it was found in samples collected in 1986 [17]. 
This allele is now present at high frequencies in An. 
gambiae s.s from Kenya, Uganda and Burundi [23-25] 
and at much lower frequencies in the same species in 
Central Africa [26,27]. 

In An. arabiensis the frequency of both the 1014S and 
1014F alleles is generally considerably lower than in An. 
gambiae [23,27,28]. In Burkina Faso, both kdr alleles 
were present in An. arabiensis although at a very low fre- 
quency (highest allele frequencies observed were 0.04 for 
1014F in 2010 in Kuinima, 0.12 for 1014S in Koupela in 
the same year). Preliminary analysis of intron data (J. 
Pinto, unpublished data) suggest that the 1014S allele in 
An. arabiensis in Burkina Faso has an independent origin 
and has not arisen by introgression from An. gambiae, 
contrary to the prediction in Kenya [23]. The source of 
selection pressure for 1014S and 1014F may differ and 
perhaps the former allele is more closely associated with 
DDT selection than pyrethroids [17]. The sample size in 
the current study was insufficient to establish phenotype: 



Badolo et al. Malaria Journal 201 2, 1 1 :232 
http://www.malariajournal.eom/content/1 1/1/232 



Page 10 of 11 



genotype correlations for the 1014S allele. However, re- 
cent functional characterization of the alternative sodium 
channel alleles supports the importance of DDT in 
selecting for 1014S by demonstrating that the L1014S 
substitution had a much greater impact on sensitivity to 
DDT than to pyrethroids [29]. 

The very high probability of surviving permethrin or 
DDT exposure in An. gambiae s.s. in Burkina Faso is at 
least partially explained by the high frequency of 1014F 
in this species. However, undoubtedly other resistance 
mechanisms are contributing towards the resistance 
phenotype. Indeed, pre-exposure to the synergist PBO 
for one hour increased permethrin mortality to >90% in 
all sites when tested in 2010 (data not shown). This is 
strongly indicative of metabolic resistance mechanisms 
and work is ongoing to investigate this further. 

The expanding geographical range and prevalence of 
pyrethroid resistance in Burkina Faso may impact on the 
efficacy of malaria control activities. There is already lim- 
ited evidence that the level of resistance to pyrethroids in 
Soumousso is impacting on the efficacy of LLINS in use in 
the field [30]. In light of this threat, the NMCP in Burkina 
Faso has elected to use bendiocarb for the IRS campaign, 
which began two years ago. Unfortunately, the ace-l R al- 
lele, which is associated with resistance to this carbamate 
insecticide, is already present in An gambiae s.s and 
appears to be spreading rapidly throughout the country. 
Although the bio-assay data do not indicate high frequen- 
cies of bendiocarb resistance at present, this situation 
should be closely monitored. 
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